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We study the possibility that the scalar mesons exist as four-quark states. The energy shift of 
two pseudoscalar mesons as a function of spatial lattice size makes a distinction between bound 
states and scattering states of four-quark states. We calculate the four-quark state in the quenched 
approximation, ignoring the two-quark annihilation diagrams and the vacuum channels. We per- 
form a calculation of pseudoscalar meson scattering amplitudes, using Nf — 2 Wilson fermion 
and plaquette/Iwasaki gauge actions. We obtain the indication that the four-quark states in the 
case of the isospin zero (7 = 0) and two (/ — 2) channels are no bound states. And we find that 
the bound energy depends strongly on pion mass rather than the ratio of pion mass to rho meson 
mass. 
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1. Introduction 

Since many light scalar mesons such as a (600), K"(800),/o(980) and ao(980) were found in 
experiments a lot of theoretical studies have been devoted to investigation of their states. For 
the structure of light salar mesons, in addition to the conventional two quark state from the quark 
model, several possibilities are proposed [Qj; four-quark states, molecular states and scattering 
states. Because the sigma meson is considered as a chiral partner of the pion in the mechanism 
of hadron mass generation, it would be interesting to investigate a roll of four-quark states in the 
mechanism. The study of four-quark states in light scalar mesons gives us insight of important 
QCD feature. 

Since Alford and Jaffe showed the possibility that scalar mesons exist as four-quark states on 
the lattice |gp, tetraquark search on the lattice started actively. The pioneer work of the sigma meson 
which is one of candidates of four-quark states was done with full QCD by SCALAR collaboration 
[Q]. They found that disconnected diagrams which contain effectively four-quark states, glueballs 
and so on make the sigma meson lighter. Recently using tetraquark interpolators, not only ground 
states but also resonance states of scalar mesons on the lattice were reported 

Following the procedure which was proposed by Alford and Jaffe [Sj, we explore the existence 
of four-quark states in the scalar channel on larger lattice with a finer lattice spacing. We also 
investigate the dependence of bound energies for four-quark states on the ratio of pion mass and 
rho meson mass (m n /m p ) and on the pion mass (m K ). 

2. Four-quark states from Lattice QCD 

We calculate four-quark correlators in two-flavor (Nf = 2) lattice QCD under the assumption 
that four-quark states of scalar mesons are consisted of bound states of the two pions. We employ 
the procedure which was proposed by Alford and Jaffe [Eh. The four-quark correlators of isospin 
zero (7 = 0) and two (7 = 2) channels are given by, 

J I=Q (t) = D{t) + hc{t) - 3A(t) + 3 -G{t) , (2. 1) 

J I=2 (t) = D(t)-C(t) , (2.2) 

where D(t), C(t), A(t) and G(t) are corresponding to the diagrams in Fig. [T[ To evaluate the four- 
quark states clearly, we carry out the calculation in the quenched approximation where two-quark, 
multi-quark and glueball states do not mix with four-quark states as intermediate states. We drop 
the contribution of the diagrams for A(t) and G(t), ignoring the two-quark annihilation diagrams 
and the vacuum channels. To obtain the bound energy 8Ej for the four-quark states, we construct 




Figure 1: The diagrams for four-quark correlators. 
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Figure 2: The physical lattice size L(= A^a) dependence of the energy shifts 8Ej = 2- In left figure (a) (right 
figure (b)), the values of m K /m p (m K ) are fixed. In both figures the lines are proportional to L~ 3 . In left 
figure (a), the three symbols are on the different L~ 3 lines, though they have the same values of m^/nip. In 
right figure (b), the two symbols which have the same values of m n are almost on the same L~ 3 line, though 
they have the different values of m K /m p . 



the ratio of the four-quark correlator Jj(t) and the pion correlator P(t) and fit it to an exponential at 
large t, 



Rr(t) 



-^exp(-S£/ ?) + ••■ 



(2.3) 



We discuss the possibility that the four-quark states exist as bound states from the Nl dependence 
of the bound energies. 

If a four-quark state is a bound state, 8Ej is negative and would approach to a negative constant 
in the large Nl region. On the other hand if a four-quark state is a scattering state, it is expected 
that 8Ej obeys the scattering formula of a two-particle state in a cube box of size L with periodic 
boundary condition in Ref. 



8E[ = Ej — 2m n 



L 3 



1+2.8373 1 ^?l+6.3752 fn * Tl 



AkL 



4kL 



+ 



(2.4) 



where Ej is the total energy, and 7} is the scattering amplitude which can be written by the scattering 
length a i as 



471(3/ 



(2.5) 



From Eq.(2.4), we find that 8Ej is proportional to L 3 in the region where the physical spatial 
lattice size L = Nlo (a is the lattice spacing) is enough large. 
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Figure 3: The physical lattice size L(= N^a) dependence of the energy shifts 8Ej = q. In left figure (a) (right 
figure (b)), the values of m K /m p (m K ) are fixed. In both figures the lines are proportional to L~ 3 . In left 
figure (a), the three symbols are on the different L~ 3 lines, though they have the same values of m^/nip. In 
right figure (b), the two symbols which have the same values of m n are almost on the same L~ 3 line, though 
they have the different values of m K /m p . 



3. Results 

In the work by Alford et al. [[J, they observed that the 1 = 2 channel is clearly the scattering 
state but found the possibility that the 1 = channel is the four-quark bound state. However to get 
the conclusive result for the four-quark bound state in the 1 = channel, further calculations on a 
larger lattice with finer lattice spacing are required. 

We calculate the four-quark correlators on a larger lattice with Nf = 2 Wilson fermion using 
configurations which are produced with both plaquette and Iwasaki gauge actions. The lattice pa- 
rameters are shown in Tables [j] and |[ We impose the Dirichlet boundary condition in the temporal 
direction on the quark fields. 

The results for / = 2 channel are shown in Fig. ||[ In Fig. § (a), the quark masses are set as 
the values of the ratio of pion mass to rho meson mass, m n /m p ~ 0.74 become close to those in 
Ref. [U (See Table @). In spite of the same values of m n /m p , the values of m K are different among 
them, which comes from the lattice artifact due to the large lattice spacings. In Fig. ^|(b), the bound 
energies are obtained under the same values of m n ~ 370 MeV. We can see that in all cases of Fig.|| 
(a) and (b) the symbols are on the L~ 3 lines, which is the same as that in Ref. [Sj. On the other 
hand the results for / = channel are shown in Fig. 0. In Fig. |3| (a), the quark masses are set as 
the values of the ratio of m^/nip become close to those in Ref. [^] (See Table |]). In Fig. |3| (b), 
the bound energies are obtained under the same values of m n ~ 370 MeV. Again we find that in all 
cases of Figs.|| (a) and (b) the symbols are on the L~ 3 lines, which is contradicted with the result of 
Ref. [§]. We can not observe the suggestion of the existence of the bound state in the 1 = channel. 

In Fig. ^ (a) the three symbols which have the same values of m K jmp are on the different L~ 3 
lines, but in Fig. [3] (b) the two symbols which have the same values of m n are almost on the same 
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Figure 4: The physical lattice size L(= N^a) dependence of the energy shifts 8Ej = q. The lines are propor- 
tional to L~ 3 . The color difference corresponds to m n difference. The green, red and blue solid lines stand 
for m% = 940, 840 and 790 MeV, respectively. The black dashed line was drawn in Fig. 5 of Ref.[^J. 



L~ 3 line. It implies that 8Ej depends strongly on m % rather than m K /m p , which is deduced from 
Eq.(15). Since the scattering amplitude is proportional to the inversion of m K , the absolute value 
of 8Ej decreases with the increasing amount of m n . 

In Ref. ^ they insisted that they found the deviation from the L~ 3 line (the black dashed line 
in Fig. |]) in the physical lattice size dependence of energy shifts in the 1 = channel. However 
we can give another explanation to it. Because of the m n differences among data plots, we can 
draw the three different L~ 3 lines (the green, red and blue solid lines) instead of one L~ 3 line (the 
black dashed line) in Fig. |j. It does not suggest that the symbols deviate from the L~ 3 line (the 
black dashed line), but suggests that the symbols are on the three different L~ 3 lines, respectively. 
It means that there are no bound states in the 1 = channel. 



4. Summary 

We have investigated a four-quark state as a bound state from the spatial lattice size dependence 
of the bound energies for a four-quark state. We found that the four-quark states in the cases of both 
1 = and 1 = 2 channels exist as no bound states, i.e., all bound energies for the four-quark states 
as a function of spatial lattice size L are on the L~ 3 lines, which is contradicted with the previous 
work [||]. The symbols which have the same values of m n /mp are on the different L~ 3 lines, on the 
other hand the symbols which have the same values of m n are almost on the same L~ 3 line. This 
results suggest that the absolute value of bound energy depends directly on m n rather than m n /m p . 

This work is a staring point of study of four-quark states in the light scalar mesons. To reach 
the conclusive results for the light scalar meson states, further investigation is needed; full QCD 
calculation with light quark masses on larger lattice size, application of four-quark interpolators, 
state-of-the-art action and so on. 
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Table 1: Lattice parameters, k is the hopping parameter. We adjust parameters, /3 and K to fix the values of 
m n /m p ~ 0.74. See Table [§ 



Gauge action 


P 


a [fm] 


K 


N[xN T 


# Conf . 


Iwasaki 


2.100 


0.278(3) [0] 


0.1666 


12 3 x 20 


150 










16 3 x 20 


120 










18 3 x20 


78 










20 3 x 20 


120 










24 3 x 20 


78 


Iwasaki 


2.300 


0.172(4) [0] 


0.1565 


12 3 x 20 


120 










16 3 x 20 


60 










20 3 x 20 


60 










24 3 x 20 


54 


Plaquette 


5.700 


0.140(4) [§] 


0.1640 


8 3 x20 


1008 










12 3 x 20 


240 










16 3 x 20 


102 










18 3 x20 


210 










20 3 x 20 


192 










24 3 x 20 


300 










36 3 x 20 


66 
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Table 2: Lattice parameters. We adjust parameters, j3 and K to fix the values of m K ~ 370 MeV. See Table 
I 



Gauge action 


j3 


a [fm] 


K 


N[xN T 


# Conf . 


Iwasaki 


2.100 


0.278(3) [§ 


0.1690 


12 3 x 20 


150 










16 3 x 20 


120 










20 3 x 20 


120 










24 3 x 20 


60 


Iwasaki 


2.300 


0.172(4) [0] 


0.1594 


12 3 x 20 


300 










16 3 x 20 


180 










20 3 x 20 


180 










24 3 x 20 


60 



Table 3: Lattice parameters in Ref. [g]. They adjust parameters, j3 and K to fix the values of m^/nip ~ 0.74. 
See Table @ 



Group 


Gauge action 


j3 


a [fm] 


K 


N[ x Nt 


# Conf. 


Gupta et al. [g] 


Plaquette 


6.000 


0.0762(8) [0] 


0.1540 


16 3 x 80 


35 


Fukugita et al. [3] 


Plaquette 


5.700 


0.162(6) [[T^] 


0.1640 


12 3 x 20 


70 


Alford et al. 


Liischer & Weisz 


1.719 


0.249(5) 




6 3 x40 
8 3 x40 




Alford et al. 


Liischer & Weisz 


1.157 


0.400(4) 




5 3 x40 
6 3 x40 
8 3 x40 





Table 4: The values of m n /m p and m K . In the upper table results of Ref. [Q] are shown. In the lower table 
results of our calculation are shown. 



data 


m n /m p 


m n [MeV] 


Gupta 




940(30) 


Fukugita 


0.740(8) 


620(90) 


Alford (a = 0.25 [fm]) 


0.756(5) 


840(11) 


Alford {a = 0.40 [fm]) 


0.756(4) 


790(6) 


Iwasaki (j3 =2.1, K = 0.1666) 


0.736(3) 


451(6) 


Iwasaki (/3 = 2.3, fc = 0.1565) 


0.746(6) 


593(16) 


Plaquette (/3 = 5.7, K = 0.1640) 


0.748(4) 


754(12) 


Iwasaki (/3 = 2.1, k = 0.1690) 


0.646(4) 


368(5) 


Iwasaki (/3 = 2.3, fc = 0.1594) 


0.545(7) 


374(11) 
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